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Allowed slepton intergenerational mixing in light of light element abundances 
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We studied allowed region on the intergenerational mixing parameters of sleptons from a viewpoint 
of big-bang nucleosynthesis in a slepton-neutralino coannihilation scenario. In this scenario, 7 Li and 
6 Li problems can be solved by considering exotic reactions caused by bound-state effects with a long- 
lived slepton. Light element abundances are calculated as functions of the relic density and lifetime of 
the slepton which considerably depend on the intergenerational mixing parameters. Compared with 
observational light element abundances, we obtain allowed regions on the intergenerational mixing. 
Ratio of selectron component to stau component, c e , is allowed in 2 x 10" 11 < c e < 2 x 10 -9 
with solving both the 7 Li and 6 Li problems. Similarly, the ratio for smuon, c M , is allowed in 
c„<5x 1(T 5 for mass difference between slepton and neutralino, which is smaller than muon mass, 
and c M < 2 x 10" 10 for the mass difference in range between muon mass and 125 MeV. We also 
discuss collider signatures of the slepton decays. We find characteristic double peaks in momentum 
distribution of event number of the slepton decays with allowed mixing parameters. Discoveries of 
the double peaks at future collider experiments should confirm our scenario. 



I. INTRODUCTION 

Long-lived Charged Massive Particles (CHAMP) are 
predicted in new physics beyond the standard model [l|. 
Such Long-lived CHAMPs induce nonstandard nuclear 
reactions in big-bang nucleosynthesis (BBN), and dras- 
tically change light element abundances [2T-l21|. By 
comparing the theoretical predictions with observational 
data, we can obtain constraints on model parameters rel- 
evant with the long-lived CHAMPs. Because the BBN is 
sensitive to lifetime from 10~ 2 sec to 10 12 sec, it should 
be one of the best tools to check an existence of the long- 
lived CHAMPs. 1 

The Minimal Supcrsymmctric Standard Model 
(MSSM) with R-parity conservation is one of the leading 
candidate providing long-lived CHAMPs. Several MSSM 
scenarios predict a long-lived stau as the next lightest 
supersymmetric particle (NLSP). An example is the 
scenario that the lightest supersymmetric particle (LSP) 
is gravitino. In this scenario, the interaction of the grav- 
itino LSP with the stau NLSP is suppressed by Planck 
scale, and hence the stau has longevity [3-[3, l30H32l |. 
Another example is an axino LSP scenario, in which 
the decay rate of the stau NLSP is suppressed because 
of loop processes or a decay constant scale F a [33U35T ] . 
Among such scenarios, the most attractive one is a 
Bino-like neutralino LSP scenario. The longevity of the 
stau NLSP is brought by tight mass degeneracy of the 
stau and the neutralino LSP. 



A remarkable feature of this mass-degenerate scenario 
is to be free from the 7 Li problem (3(|. The 7 Li problem 
is a discrepancy between the observed abundance of 7 Li, 
e.g., Log 10 ( 7 Li/H) = -9.63 ± 0.06 [37|, and a theoreti- 
cal one predicted in the standard BBN, Log 10 ( 7 Li/H) = 
—9.35 ± 0.06 [2l[. (There is also a severer observational 
limit which worsens the fitting [HI].) 2 In order to explain 
the observed abundance of the dark matter by the Bino- 
like neutralino LSP based on a thermal relic scenario, 
the so-called stau coannihilation mechanism is required 
to work well (45j . This mechanism requires of order or 
smaller than 1% degeneracy in mass between the neu- 
tralino LSP and the stau NLSP. Due to this tight de- 
generacy, the stau can be long-lived j4(| |47j and trigger 
exotic nuclear reactions at the BBN era. 3 One of such ex- 
otic reactions is an internal conversion process in a bound 
state of the stau and 7 Li ( 7 Be) nuclei [ToL [TT|. In the in- 
ternal conversion process, 7 Li ( 7 Be) can be destroyed into 
lighter nuclei and their abundances are fitted to the ob- 
servational ones. Therefore, in the mass-degenerate sce- 
nario, the observed abundances of the dark matter and 
all light elements are simultaneously realized. 

It is further worth examining the BBN in the mass- 
degenerate scenario. In general, the MSSM has sources 



1 See also Refs. I22I429H for long-lived neutral particles. 



This discrepancy can not be solved even for corrections of cor- 
responding cross sections of nuclear reaction [39l , l4fj . Even if 
we consider nonstandard astrophysical models including diffu- 
sion process es |41| . |4~3 ]. it would be difficult to fit all of data 
consistently |43l . See also [43 for a recent work related to light 
element nucleosynthesis in accretion flows. 

See also [48H52II for another cosmological constraints on 
CHAMPs. 
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of intcrgcncrational mixing, though wc had assumed the 
lepton flavor conservation in our previous works. With 
the intergenerational mixing, the mass eigenstates are 
linear combinations of flavor eigenstates. As we will 
discuss later in detail, the intergenerational mixing re- 
duces the number density of the slepton via two types 
of processes. Success and failure of solving the 7 Li prob- 
lem strongly depends on the number density of the slep- 
ton. Furthermore, major and minor of modifications on 
the light element abundances are also determined by the 
number density of the slepton. Thus, the modified abun- 
dances of the light elements can be predicted by the num- 
ber density that is a function of mixing parameters. As 
a result, we can constrain the mixing parameters from 
the consistency between the modified and the observed 
abundances of the light elements. 

The aim of this work is to predict the slepton mixing 
parameters in the mass-degenerate scenario. In this sce- 
nario, the exotic nuclear reactions catalyzed by the long- 
lived slepton induce both a fusion process Q and destruc- 
tion processes through internal conversions for 7 Li and 
7 Be [lMlHl, and spallations for 4 He These pro- 
cesses can modify the abundance of the light elements. In 
fact, 6 Li is produced via the former process after forming 
a bound state of the slepton with 4 He. A fit b y th e ob- 
servational abundance 6 Li/ 7 Li = 0.046 ±0.022 [54| gives 
an allowed region in a parameter space of the mixing pa- 
rameters. In addition, both deuterium (D) and tritium 
(T) (or 3 He after its decay) are produced nonthcrmally 
by the latter processes. By adopting recent observational 
constraints LVH = (2.80±0.20) x 10~ 5 [H[ and 3 Hc/D < 
0.87 + 0.27 [56j, an upper bound on the number density 
of the slepton is obtained as a function of its lifetime. 
Thus, the upper bound on the number density results 
in a lower bounds on the mixing parameters of the slep- 
ton. On the other hand, a sufficient number density of 
the slepton at the bound state formation with 7 Li ( 7 Be) 
is required for solving the 7 Li problem. Then, an upper 
bound on the mixing parameters is obtained. Hence, in- 
triguingly, the mixing parameters of the slepton are pre- 
dicted with pinpoint accuracy from both above and below 
in light of the current observed light element abundances. 

This long-lived slepton scenario can be also exam- 
ined at terrestrial experiments. For example, due to 
the longevity longer than "The first three minutes of 
the universe", a small number of the long-lived slepton 
would decay in detectors 4 . Some of them decay into 
the neutralino LSP and electron (or muon) via the in- 
tergenerational mixing. Combined with the predictions 



on the mixing parameters using the light element abun- 
dances, we find monochromatic and diffuse spectrum in 
momentum distributions of the final electron (or muon). 
These spectrum are characteristic signatures of the mass- 
degenerate scenarios, and therefore we can distinguish 
our scenario from others. 

This paper is organized as follows. In the next section, 
we discuss the relation between the number density of 
the long-lived slepton and the intergenerational mixing 
in the mass-degenerate scenarios. Then, we present a set 
of Boltzmann equations for calculating the number den- 
sity as a function of the mixing parameters. In Sec. IIII1 
we analytically estimate bounds on the mixing parame- 
ters. Then, in Sec. IIV1 we show our numerical results. 
Allowed region on the intcrgcncrational mixing parame- 
ter is shown from the viewpoint of observed light element 
abundances. In Sec. [Vj applying the results, we discuss 
collider implications of this long-lived slepton. Finally, 
we summarize our results and discuss the prospects. 



II. INTERGENERATIONAL MIXING AND 
NUMBER DENSITY OF LONG-LIVED SLEPTON 

Long-lived slepton leads some kinds of exotic nuclear 
reactions and changes light element abundances. Pre- 
dicted light clement abundances depend on the number 
density of the slepton at the BBN epoch. Interestingly, 
the number density is quite sensitive to the intergener- 
ational mixing parameters of the slepton in the mass- 
degenerate scenarios. Therefore, in order to constrain 
the intergenerational mixing, the number density has to 
be accurately calculated as a function of the intergener- 
ational parameters. 

A calculation of the number density of the slepton con- 
sists of two steps. In the first step, the total number 
density of SUSY particles is calculated at a time of their 
chemical decoupling. After the chemical decoupling, the 
number density of the slepton continues to evolve till the 
temperature T ~ 5m via scattering off the cosmic ther- 
mal background. Here dm is the mass difference between 
the slepton NLSP and the lightest neutralino LSP. Fur- 
thermore the number of the slepton is reduced by its own 
natural decay Then, as the second step, the evolution of 
the number density is acquired by solving relevant Boltz- 
mann equations implementing the scatterings and the de- 
cays. We explain the calculation of the number density, 
and then present a set of Boltzmann equations for the 
calculation. 



4 Collider signatures of the long-lived slepton in the neutralino 
LSP scenario is studied in case that its lifetime is shorter than 
the beginning of BBN l57ll58l in collider experiments. Measuring 
the lifetime and the mass difference between the slepton and the 
neutralino LSP, the intergenerational mixing will be determined. 
Combining with the other works on collider signatures of long- 
lived, the models can be discriminated. 



A. First step: Total number density of SUSY 
particles 

As the first step, we calculate the total number density 
of SUSY particles. The total number density is equal to 
the relic number density of the neutralino dark matter 
in the current universe because all SUSY particles decay 
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FIG. 1: Lifetime of the slepton as a function of 8m for mj = 
350GeV and sin<9 e = 0.8. Solid line is the lifetime of the 
pure stau, and others are the lifetime of slepton with each 
intergenerational mixing parameters. 



into the neutralino LSP in the end under the R-parity 
conservation. Calculation of the relic density of the neu- 
tralino dark matter in the stau-neutralino coannihilation 
scenario has been developed so far [III |59|. Assuming 
the intergenerational mixing is not so large and the main 
component of the slepton is stau, the calculation method 
is applicable for our calculation. That is, cross sections of 
coannihilation processes are firstly calculated, and then 
we solve the Boltzmann equations for sum of the number 
densities of SUSY particles. 

After the chemical decoupling, the total number den- 
sity of SUSYparticles remains the value at the freeze out. 
This value does not have correlation with the intergen- 
erational mixing. This value is required as the initial 
condition for calculating the slepton number density as 
a function of the intergenerational mixing. 



B. Second step: Ratio of the number density 
between the slepton to the neutralino 

A key ingredient for constraining the mixing parame- 
ters is the evolution of the number density of the slepton 
after the chemical decoupling of SUSY particles. 

Thermally-produced slepton decays according to its 
lifetime after the chemical decoupling. Hereafter we as- 
sume that the mass difference 8m is smaller than the mass 
of tau lepton. In this situation, the decay of slepton into 
the neutralino LSP and a tau lepton are forbidden kine- 
matically. In this case the slepton becomes long-lived 

M- 

The decay of the slepton is categorized into two types. 
The first one is 3- or 4-body final state processes through 
the stau component; 



I ->Xi+Vt+1 +"l (Zee,//). 



(1) 



In the mass-degenerate scenario, the decay rates of these 
modes considerably depend on 8m (see Fig. Q] on this 
work and Fig. 2 in Rcf. [47|). The second one is 2-body 
final state processes through the selectron or smuon com- 
ponent; 



(2) 



The 2-body decay process are dominant when the mixing 
parameters are not very small so that the decay rates 
are larger than those of 3- and 4-body decays. On the 
contrary, when the intergenerational mixing is tiny, 3- 
or 4-body decay processes dominate over 2-body decay 
processes in spite of a tight kinematical suppression. The 
constraint on the mixing from these decay processes must 
be therefore studied carefully as both functions of 8m and 
the mixing parameters. 

Exchange processes by scattering off the thermal back- 
ground also play an important role in constraining on the 
mixing parameters. Even after the chemical decoupling 
of SUSY particles, the ratio of the number density of the 
slepton to the neutralino continues to evolve via exchange 



processes 



^+7^X1+/^, 



(3) 



Note that the process Z7 O Xi e must not be included. 
The decay mode I x?e is always open in the presence 
of the intergenerational mixing, and the process I"/ 
X±e should be counted as one of the corrective part for 
this decay (or inverse decay) mode with soft emission (or 
absorption). Thus, including the process Z7 O Xi e m 
addition to this decay mode results in a double counting. 
Similarly, when the mode I — > XiM i s open, then the 
corresponding exchange process Z7 -O- x?M also must be 
taken out. 

These exchange processes keep the slepton and the 
neutralino in kinematical equilibrium. Due to this kine- 
matical equilibrium, the number density ratio between 
the slepton and the neutralino maintain the Boltzmann 
distributions. The number density of the slepton for a 
fixed dm therefore decreases as the temperature of the 
Universe decreases; 



-N = e- 5m ' T 



N 



2(i + e - Sm / T ) ' 



(4) 



where T is the temperature, and N is a sum of the densi- 
ties of the slepton and the neutralino that are obtained in 
the first step (see Sec. Ill A| ). The kinematical equilibrium 
is broken when the Hubble expansion rate overwhelms 



5 Although other exchange processes exist, those contributions are 
negligible due to weak interactions. 
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the reaction rate of the exchange processes. Hence it 
is important to know when the Hubble expansion rate 
overwhelms the reaction rate. We numerically solve a 
set of Boltzmann equations to obtain accurate resultant 
number density of the slepton. 



light element abundances. Comparing them with obser- 
vational values, we can constrain or predict the mixing 
parameters. 



C. A set of Boltzmann equations for relic number 
density of the slepton 

Now we present a set of the Boltzmann equations to 
calculate the number density of the slepton as a function 
of the intergencrational mixing parameters. The Boltz- 
mann equations describing evolutions of the number den- 
sities of the negatively charged slepton , the positively 
charged slepton nr + , and the neutralino are, 



dn i- 
dt 



i,X,Y 



- E [( r W.. fi i-"( r )i- fi x4V- 



X,Y.. 



dnj. 



(5) 



l+X' H ' l Y 



i,X,Y 



- E [( T )xXY...n I+ - (r) r+ n x n^ny 9 .. 



X,Y,... 



(0) 



dn 



dt x 



E \{ Cjlv )iY n X n 'x ~ {°' v )xX n i n 
i,X,Y 

E [( T )i n X n X n Y - - { r )xXY..J 



,X.Y,. 



(7) 



where t is time and H is the Hubble parameter. In the 
right hand side, (cr'v) is the thermal averaged cross sec- 
tions of the exchange process in Eq. ([3]), and (r) is the 
thermal averaged decay rate (or inverse decay rate) . The 
index i represents relevant SUSY particles, and indices 
X and Y represent the SM particles, respectively. Sub- 
scripts to both (cr'v) and (r) denote the final states of 
each process. In the equations, we assumed that the rel- 
evant SM particles are thermalized 6 . 

In the presence of the intergencrational mixing, (cr'v) 
and (r) depend on the mixing parameters. Solving a set 
of the Boltzmann equations above, the number density of 
the slepton at the BBN epoch is obtained as a function 
of the mixing parameters. By using obtained number 
density, we perform the BBN calculations and predict 



III. ANALYTICAL CONSTRAINT ON THE 
INTERGENERATIONAL MIXING 



Before showing our numerical results, we give analyti- 
cal estimations for bounds on the mixing parameters. We 
expand the lightest slepton in terms of the flavor eigen- 
states, e, jl and f without the intergenerational mixing, 



f = cos 8 f f L + sin f f R , 



(8) 



where fi, and fa are left-handed and right-handed slep- 
tons, and Of are the mixing angles of these sleptons, re- 
spectively. The coefficients c e , c M , and c T are the mixing 
parameters. These parameters are in ranges from zero to 



unity and satisfy the relation c 2 



= 1. In this 



work, we assume c e 



and c T to be real and positive 



parameters for simplicity. 

These mixing parameters are constrained from above 
by two requirements. The first requirement is the suffi- 
cient longevity of the slepton, and the second one is the 
sufficient number density at the moment of decoupling 
from the exchange processes. For the first requirement, 
the lifetime of the slepton must be longer than O(10 3 s) 
to solve the 7 Li problem. The slepton with this lifetime 
can form a bound state with 7 Li and 7 Be and des troy 
these nuclei via the internal conversion processes [TJl fllj . 
The lifetime is depicted as a. function of 8m in Fig. [1] for 
to,- = 350GeV and sin6> e = sin 6^ = 0.8. A curve labeled 
"stau" is the lifetime for a pure stau, and other curves 
arc the lifetime for the slepton with the mixing parame- 
ters labeled near the curves. The lifetime varies by a few 
order of magnitudes as 8m and c e ^ vary. Due to the kinc- 
matical thresholds for the decays, dependence on both c e 
and Cfj, is different for 8m > m M or 8m < m^. There- 
fore, the lifetime should be studied in each 8m region to 
obtain bounds on the mixing parameters. 

To obtain the bounds analytically, we consider two sim- 
ple cases in which either c e = or c M = for simplicity. 
For c M — case, the slepton decays into an electron and 
the lightest neutralino (2-body decay). Then the lifetime 
of the slepton is approximately given by 



Tj(l^xl + e) 



8tt 



1 



g 2 tan 2 0w (8m) 2 cos 2 6 e + 4 sin 2 



(9) 



Vability of this assumption is checked in Ref. [53 



where mj is the mass of the slepton. By inserting cos ( 



5 



0.6 for a reference value, the lifetime is expressed as 



O.lGcV 



Sm 



1.34 x 10 3 



3(0.6) 5 



III; 



300GeV 

- 10 -10.B 



4-3 cos 2 9 e 



Requiring tj > O(10 3 s), we obtain the bound on the 
selectron component, c e < 3.2 x 10~ n for Sm = 0.1 GeV 
and mj = 300 GeV. Similarly, for c e = case, the lifetime 
is approximately given as 



10 d 



io- 



0.12GcV 



Sm 



mj 


/ 25 
V 36 


'0.12GeV~ 




300GcV 


Sm 





-1/2 



(10) 



where cos^ = 0.6 is used. From Eq. (fT0|) . the upper 
bound on c M is obtained as 3.2 x 10 -11 for Sm = 0.1 GeV 
and mj = 300 GeV. 

We can derive another upper bounds on the mixing 
parameters from the second requirement that the num- 
ber density should be sufficient at the decoupling time 
from the exchange processes to destroy the 7 Li and 7 Be 
nuclei. Omitting smuon component again, the exchange 
processes are le O xll ancl ' T Xi7 through selec- 
tron and stau component. The former exchange process 
keeps working even after the freeze out of the latter pro- 
cess. This is because densities of electron and photon 
are still quite large in thermal bath for T < m T . Then, 
as explained in the previous section, the slepton number 
density continues to decrease. The reaction rate with 
the electron is proportional to c\. Hence the larger the 
mixing is, the more the number density decreases. The 
reaction rate of le <->■ Xil must be suppressed to ensure a 
sufficient number density of the slepton. The bound on 
c e is estimated by comparing this reaction rate with that 
of processes with tau. 

In the absence of the intcrgcnerational mixing or pure 
stau case, freeze out of the exchange processes (/ ± 7 O 
X°t ± and x?7 «-> ^ ± t =f ) occurs at T ~ 70MeV, and 
is almost independent of both Sm and mj (53[. Then 
the number density of the slepton manages to be a suf- 
ficient amount. This fact suggests that the reaction rate 
of le O XiY must be smaller than that of the pro- 
cesses in the absence of the intergcncrational mixing at 
T = 70McV. Parameterizing the cross section of this pro- 
cess as (a'v) e = c 2 e (a'v) T , the ratio of these reaction rates 
at T = 70MeV is given by 



(a'v) e Y i Y e 
(a'v) T Y i Y T 



(1.08 x 10 9 )c 2 



(11) 



at Yi and Y e ^ T are the yield value of the slepton and 
electron/tau, respectively. The mixing parameter c e is 
therefore required to be < 3 x 10 -5 . For c e = case, 
the reaction rate of the process l\i -O- x?7 is compared 
with that without the mixing. Parameterizing the cross 



section of the process as (<r'v)^ — c 2 l (a'v) T , the ratio of 
the reaction rates is given by 



(<j'v)uYjY u , 7N 9 

Ut# * (9 ' 93 x 10 )c - 



(12) 



at temperature T = 70MeV. From Eq. ([T2]). c M must to 
be smaller than 10 -4 . 

As we will see later, we can constrain c e from below by 
the relic abundance of light elements. Thus we will get 
an allowed region for c e . 



IV. NUMERICAL RESULTS 

We are now in position to numerically search for an 
allowed region of the slepton mixing parameters. First 
we see that the upper bounds on c e and c M obtained in 
the previous section are in a good agreement with nu- 
merical analysis. Then we numerically compute reaction 
networks of light elements including the exotic nuclear re- 
actions, and find parameter regions as a function of Sm 
allowed by the observational light element abundances. 



A. Constraint on slepton intergenerational mixing 



First we show bounds on the mixing parameter c M for 
Sm < rria- Based on previous works we fix Yr = 
2 x 10~ 13 for the sufficient number density to solve the 
7 Li problem. 

Figure [2] shows the temperature evolution of the yield 
values of the slepton for each values of c^. The value 
of c M is shown near corresponding curves. From top 
to bottom panels, the mass differences are taken to be 
Sm = 20 MeV, 60 MeV, and 100 MeV, respectively. The 
curves tagged by "stau" represent the yield values with- 
out the mixing, and the curves tagged by "Equilibrium" 
are the yield values of the slepton in kinetic equilibrium. 
In shaded region, number density of the slepton is insuf- 
ficient to solve the 7 Li problem. In each panel, we took 
c e = 0, because bound on c e is c e < 10 -10 and hence 
contributions of c e are negligible when the exchange pro- 
cesses are still working. As is expected, larger value of 
keeps yield values of the slepton to be in kinetic equilib- 
rium and leads smaller number densities of the slepton. 

The bound on c M < 5 x 10~ 5 for Sm = 20 MeV 
(Fig. [3](a) ) accurately reproduces our estimation in the 
previous section. For larger Sm, the bounds on c M be- 
come more stringent; c M < 2 x 10 -6 for Sm = 60 MeV 
(Fig. Mb)), and c M < 2 x 10~ 7 for Sm = 100 MeV 
(Fig. [2Jc)). This is understood as follows. The cross sec- 
tions for the processes ([3]) are not so dependent on Sm 
for Sm < m M . They are more dependent on the muon 
number density and hence the temperature. It means 
that the decoupling temperature of these processes are 
almost same as long as is same. Then to have a suf- 
ficient yield, these processes must decouple earlier for 
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Fig.[3]shows bounds both on c p and c e for 5m = 106MeV, 
114MeV, 122MeV, respectively. Values near each curves 
are the lifetime of the slepton. In shaded region, the 
lifetime of the slepton is shorter than 1000s, and hence 
the parameters in this region are excluded. In left-side 
and down-side region in each panel, 4-body decay process 
is dominant because the mixing parameters are so small. 
This result is consistent with the result in Fig. [TJ The 
upper bounds on c M obtained here are c M < 2 x 10~ 10 for 
5m = 106MeV, which is just above the threshold of the 
decay of the slepton into a muon, and c M < 6 x 10~ n for 
5m = 114MeV, c„ < 2 x 10" 11 for 5m = 122MeV. 
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FIG. 2: Evolution of slepton yield values Yj for each smuon 
component c M , which values are attached on corresponding 
curve. Curves tagged by "stau" are yield values of pure stau, 
and curves tagged by "Equilibrium" are yield values of the 
slepton in kinetic equilibrium. In shaded region, number den- 
sity of the slepton is insufficient for solving the 7 Li problem. 



larger mass difference. See cq.(|l]). Thus less intergenara- 
tional mixing is allowed for larger mass difference. 

In the case of 5m > m^, more stringent bounds on 
c M are derived from the requirement on lifetimes of the 
slepton. To solve the 7 Li problem, the bound states of 
the slepton and 7 Be ( 7 Li) have to be sufficiently formed. 
This formation requires the lifetime of the slepton to be 
at least 1000s, preferably 2000s (see Fig. 1 in Ref. [3). 



B. Allowed regions by big-bang nucleosynthesis 

By numerically solving the Boltzmann equations, we 
obtain time evolutions of the bound states such as ( 4 He 
( 7 Li l~), and ( 7 Be l~) including charge-exchange re- 
actions [l6|, [5l[. Once those bound states are formed, 
the elements are immediately destroyed through inter- 
nal conversion processes for 7 Li and 7 Be, and spallation 
processes for 4 He, induced by bound-state effects. (lOr - 
[I2], [5^|. Then by solving coupled equations of the re- 
actions including those nonstandard processes, we can 
obtain final abundances of the light elements. We adopt 
a value of baryon to photon ratio 77 = (6.225 ± 0.170) x 
10" 10 (68% C.L.) which was reported by the WMAP 
satellite (rJoj . 

7 Be and 7 Li are efficiently destroyed through internal 
conversion processes ( 7 Be l~) — > 'Li + v\ + Xi with a 
following standard process 7 Li + p — ^ 4 He+ 4 He, or an- 
other subdominant internal conversion process ( 7 Li i - ) 
— > He + vi + x . Then we can compare this theoreti- 
cal value with the observational abundance. It is notable 
that most of primordial 7 Li is produced by primordial 
7 Be through its electron capture at a much later time 
for the adopted value of rj. We have also included the 
production process of 6 Li through the bound-state effect 
( 4 HeP) + D -> 6 Li + l~UM- 

On the other hand, D and 3 He are nonthermally pro- 
duced by the 4 He spallation processes ( 4 He^~) — > T +n 
+ Xi + n, ( 4 Hcf-) D + 2 n + x°i + v u ( 4 Hc[") 
p + 3 n + Y? + vi, and subsequent standard processes 
for n and p [531 ■ Note that 7 Li or 7 Be is also secondar- 
ily produced by those nonthermally-produced energetic 
T and 3 He, which may worsen the 7 Li problem partly in 
the parameter spaces. 

In Fig. 2] we plot regions allowed by observational 
light element abundances in 5m—Yj_ planes in cases of 
c e = 10~ n (top panel), 5 x 10~ u (middle panel), and 
10~ 10 (bottom panel). Lines of the constraints are plot- 
ted for D/H, 3 He/D, and 6 Li/ 7 Li at 2cr. An exception 
is 7 Li/H whose lines are denoted by both dotted lines at 
2a, and solid lines at 3cr. The theoretical curve of the 
relic abundance is also plotted as a thick solid line. For 
simplicity, here we have assumed = 0. 

In Fig. [51 we also plot allowed regions in the 5m— c e 



7 



* 5X1CT 11 
15 

= 

o 

a* ixio -11 

o 

U 

= 

o 
H 
u 

50 1X10" 12 



Insufficient longevity 




mass difference 106MeV 



(a) 



1 X 10" 11 
5X 10"" 



1 X 10" 



3X 10" 



s 

a 

o 

a* i x io" n 

o 

U 

s 
© 

H 
u 

"33 

1X 10-12 





mass difference 11 4MeV 




ioooX 


Insufficient longevity 


' lSOOs^v^^ 








(b) 



1 X 10" u 



1 X 10" 



3X 10" 



„ 5X10" 11 

= 

Hi 

ca 
o 

g 1 1 X 10 -11 

© 

U 

= 
o 
- 

u 

!/2 1X10 -12 





mass difference 122MeV 




Insufficient longevity 








(c) 



1X10" 11 1X10" 10 3X10" 10 

Smuon Component 



FIG. 3: Bounds both on c M and c e for 8m > m M . Shaded re- 
gion is disfavored in light of solving the 7 Li problem, wherein 
the slepton decays before forming a bound state with 7 Be and 
7 Li. 



plane by using the theoretical value of the relic abun- 
dance for the negatively charged slepton at each point of 
the plane. Each panel is plotted with same parameters. 

The region on the bottom-left side of thick dotted 
line is free from 7 Li problem at 3cr, and the region on 
the other side of the line is excluded due to insufficient 
density of the slepton. The regions between two blue 
solid lines are allowed for 6 Li/ 7 Li at 2a. The 4 He spal- 
lation processes exclude shadowed, light shadowed, and 
dark shadowed regions due to over-productions of D/H, 
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The 4 He spallation processes exclude shadowed, light shad- 
owed, and dark shadowed regions due to over-productions 
of D/H, 3 He/D, and 7 Li/H, respectively. The region on the 
bottom-left side of left blue solid line is excluded by 6 Li over- 
production. The white region in top panel (bottom panel) is 
the parameter space which is consistent with all of observa- 
tional light element abundances including 7 Li (both 7 Li and 
6 Li). 



3 Hc/D, and 7 Li/H at 2er, respectively. The region on 
the bottom-left side of left blue solid line is excluded by 
6 Li over-production. Thus the white region in top panel 
(bottom panel) is the parameter space which is consis- 
tent with all of observational light element abundances 
including 7 Li (both 7 Li and 6 Li). 
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FIG. 6: Event number distribution of the decay of the long- 
lived slepton as a function of momentum of final state elec- 
tron. We assumed detector resolution to be lMeV, and nor- 
malized number of 2-body decay to 10 4 events in a bin. 



V. COLLIDER SIGNATURE 

Finally, we have a short discussion on collider sig- 
natures of the long-lived slepton. The long-lived slep- 
tons stopped in a detector decay into 2-body final state 
(I — > Xi e ) an d 4-body final state (I —> XiV T ev e ). From 
the results in previous section, we can predict the ratio 
of the number of events between these decays. One of 
the way to confirm our scenario is to search for double 
peaks in electron momentum distribution in the decay of 
the long-lived slepton, and count the event number ratio. 

Figure |JJ] shows expected event distributions of the de- 
cay as a function of electron momentum. The event dis- 
tributions of both 2-body and 4-body decays are plotted 
in log scale in top panel, and bottom panel plots the dis- 
tribution of only 4-body decay in linear scale. We took 
mj = 350GeV, 5m = 0.07GeV, c e = 2 x 1CT 11 , and 
c M = 0. This parameter set is found in allowed region 
(see Fig.[S]). We assumed detector resolution to be lMeV, 
and normalized number of 2-body decay to 10 4 events in 
a bin. One can see in the top panel that there is a peak 
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at the momentum of 0.07GcV. The peak corresponds to 
2-body decay because the momentum is the same value 
of 5m. On the other hand, in the bottom panel, one 
can see that 4-body decay also shows a peak around 0.2 
GeV. This peak is broader than that of the 2-body de- 
cay. Number of event of 2-body decay is almost 100 times 
more than peak number of 4-body decay Thus, the event 
distribution shows a double peak structure, and this is a 
characteristic feature for mass-degenerate scenario. 

VI. SUMMARY 

We have considered a scenario of the MSSM where the 
slepton NLSP is long-lived due to a small mass difference 
between the NLSP and the ncutralino LSP, and studied 
the effects of the intergenerational mixing of sleptons on 
big-bang nucleosynthesis. In this scenario, the so-called 
internal conversion processes occurs in the bound states 
between the slepton NLSP and light nuclei in BBN. Then, 
the 7 Li and 6 Li problems can be solved simultaneously 
when the slepton NLSP is enough long-lived and its num- 
ber density is sufficient at the era of BBN. We have ana- 
lyzed the yield value and the lifetime of the slepton with 
the intergenerational mixing as well as the relic abun- 
dances of the light elements including the internal con- 
version processes in BBN. 

In Sec. Ill, we have calculated the yield value and the 
lifetime of the slepton NLSP and derived upper bounds 
on the mixing parameters c M . The upper bounds are 
obtained by requiring the yield value and the lifetime to 
be Yj ~ 2 x 10~ 13 and tj > 10 3 sec. to solve the 7 Li 
problem. In the case of 5m < m^, the upper bounds 
are given as c M < 5 x 10~ 5 , 2 x 10~ 6 and 2 x 10~ 7 for 
8m = 20, 60 and 100 MeV, respectively. On the other 
hand, in the case of 5m > m^, the upper bounds, c M < 
2 x 10~ 10 , 6 x 10~ n and 2 x 10~ n , are given for 5m = 
106, 114 and 122 MeV, respectively. 



We have also analyzed the the relic abundances of the 
light elements by solving the Boltzmann equations tak- 
ing into account all the exotic processes listed in SeclIVI 
We derived an allowed region on c e by adopting 2a un- 
certainties for D/H, 3 He/D and 6 Li/ 7 Li, and 2a and 3a 
uncertainty for 7 Li/H. We found that, assuming = 0, 
the 7 Li and 6 Li problems can be solved simultaneously 
in the range of 2 x 10~ n < c e < 2 x 10 -9 . These results 
are consistent with the estimations given in Sec. II. 

In the end, we have discussed the decays of the slepton 
in a case of 5m = 0.07 GeV and c e = 2 x 10 -11 as an illus- 
trating example. We found that the 2-body decay shows 
a sharp peak at the momentum of an outgoing electron 
close to while the 4-body decay shows a broad peak in 
the momentum distribution. The sharp peak is found at 
the momentum equal to 5m, and the expected number of 
events for 2-body decay is 100 times larger than that for 
4-body decay. These features as well as heavy charged 
tracks are unique in the degenerate-mass scenario. In 
the case that lifetime of the long-lived slepton is 10 3 sec, 
even present working detectors can stop the produced 
long-lived sleptons and accumulate ~0.5% of them ap- 
proximately [6l|. If exotic charged tracks of long-lived 
particle are discovered, new facility may be constructed 
with higher accumulation efficiency and high momentum 
resolution. Discovery of unique signatures of the decay 
in addition to charged tracks of long-lived particle will 
be a strong evidence of this scenario. 
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